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Low-energy dipole excitations are analyzed for the stable isotopes *''Ca and **Ca in the framework 
of the Skyrme-second random-phase approximation. The corresponding random-phase approxima- 
tion calculations provide a negligible strength distribution for both nuclei in the energy region from 5 
to 10 MeV. The inclusion and the coupling of 2 particle-2 hole configurations in the second random- 
phase approximation lead to an appreciable dipole response at low energies for the neutron-rich 
nucleus **Ca. The presence of a neutron skin in the nucleus ''^Ca would suggest the interpretation 
of the low-lying response in terms of a pygmy excitation. The composition of the excitation modes 
(content of 1 particle- 1 hole and 2 particle-2 hole configurations), their transition densities and their 
collectivity (number and coherence of the different contributions) are analyzed. This analysis indi- 
cates that, in general, these excitations cannot be clearly interpreted in terms of oscillations of the 
neutron skin against the core with the exception of the peak with the largest B{E1) value, which 
is located at 9.09 MeV. For this peak the neutron transition density dominates and the neutron 
and proton transition densities oscillate out of phase in the internal part of the nucleus leading to 
a strong mixing of isoscalar and isovector components. Therefore, this state shows some features 
usually associated to pygmy resonances. 

PACS numbers: 21.60.Jz,21.10.Re 
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I. INTRODUCTION 

The evolution of the low-lying dipole response as a 
function of the isospin asymmetry has been extensively 
analyzed both experimentally and theoretically in sev- 
eral stable and unstable nuclei. A recent review about 
the main theoretical results and experimental measure- 
ments can be found in Ref. [11]. Measurements on the 
halo nuclei ^He, *^He, "Li, ^^Be, ^^C and ^B (Refs. [l|,[2 
and references therein) have shown the presence of a sig- 
nificant dipole strength at very low energy (for instance, 
in the case of "'^^Li, below 4 MeV). In these light nuclei, 
the development of a strong dipole response at low en- 
ergies can be related to the extremely small value of the 
separation energies (the systems are very weakly bound 
and the last occupied neutron states are close to the con- 
tinuum states). These low-energy states are not collec- 
tive excitations and have mainly a single-particle char- 
acter Q. The same kind of picture (individual excita- 
tions of the weakly bound last occupied neutron states) 
is currently provided by most of the available theoretical 
models in the description of the low-lying excitations in 
neutron-rich Oxygen isotopes. Low-energy dipole exci- 
tations have been observed for the isotopes 16^22q ^^^ 
the corresponding data are reported in Refs. [^. From 
the experimental point of view, the character of these ex- 
citations (collective or single-particle) has not yet been 
clearly elucidated. 

For heavier stable and unstable nuclei, the develop- 
ment of a pygmy dipole response in neutron-rich systems 
is currently related to the formation of a thick neutron 
skin at the surface of the nucleus: the low- lying dipole 



modes are interpreted in terms of oscillations of the skin 
against the core composed by both neutrons and pro- 
tons. Experimentally, low-lying Ei states have been mea- 
sured in several medium-mass and heavy nuclei such as, 
for example, 40,44,48^^ gQ, the tin isotopes "^Sn 0, 
"6,i24sn a i30Sn and i^zsn @, 204,206,207,208pb ^^ 
and^SNi [nl. 

From the theoretical point of view, pygmy resonances 
have been analyzed with several models. Some exam- 
ples are the relativistic and non-relativistic (Q)RPA ap- 
proach (see, for instance, Ref. [l|] and references therein, 
and Refs. [i4|-[3), the particle-phonon-coupling models 
(Ref. [l| and references therein), a semiclassical coupled- 
channels approach for Sn isotopes [13, [l3| , hydro dynam- 
ical models [l9|, the phonon-damping model [20|, and 
the so-called Extended Theory of Finite Fermi Systems 
(ETFFS) [21[. Discrepancies among the different theo- 
retical predictions are found concerning in particular the 
collective character and the fragmentation of the low- 
lying modes. 

The widths and the fragmentation of the excited modes 
in a many-body system cannot be described within the 
standard random-phase approximation (RPA) which can 
only account for the so-called Landau damping (related 
to single-particle degrees of freedom). It is well known 
that, to describe widths and fragmentation, the single- 
particle degrees of freedom have to be coupled with 
more complex configurations (collective coordinates or 
multiparticlc-multihole configurations) within a beyond 
mean-field model. Among the different beyond mean- 
field models that allow one to describe, at least par- 
tially, the width and the fragmentation of the exci- 



tation modes, the second random-phase approximation 
(SRPA) is a powerful theoretical tool where the cou- 
pling with 2 particle- 2 hole (2p2h) configurations is in- 
cluded within an RPA-like formalism. In this way, the 
so-called spreading widths can be described together with 
the Landau damping. Escape widths are missing if the 
coupling to the continuum is not included. To include 
higher multiparticle-multihole configurations different di- 
rections may be followed |22h24| . 



Due to the heavy numerical effort required, the SRPA 
equations have been often solved resorting to some ap- 
proximations, namely the SRPA equations have been re- 
duced to a simpler second Tamm-Dancoff model (Le. the 
matrix B is put equal to zero, see for instance [25|-|28|) 
and/or the equations have been solved with uncorrelated 
2p2h states in the so-called diagonal approximation (29l - 
l35l |. Recently, full SRPA calculations have been per- 
formed for some O and Ca isotopes [3a, 123] ■ In particu- 
lar, in Ref. [33 calculations with the density-dependent 
Skyrme interaction have been performed adopting two 
currently used approximations for treating the rear- 
rangement terms of the residual interaction appearing in 
beyond-RPA matrix elements. The two approximations 
consist in cither neglecting these rearrangement terms or 
treating them with the standard RPA procedure. Impor- 
tant differences have been found between the correspond- 
ing two sets of results. The same authors have addressed 
this point in a more recent work [38| where a procedure 
to derive the expressions of all the rearrangement terms 
within the SRPA framework has been presented and ap- 
plied to calculations for the nucleus ^^O. In this first ap- 
plication, the importance of the proper treatment of the 
rearrangement terms in SRPA for the description of the 
fragmentation of the excited modes has been shown. 



In this work, we employ the implemented code where 
the full rearrangement terms have been included to treat 
the low-lying excitation spectrum of the stable isotopes 
""^'"^^Ca within the Skyrme-SRPA model. Medium-mass 
Ca isotopes are chosen as intermediate cases between 
light nuclei where the low-lying dipole excitations have 
mainly a single-particle character and heavier nuclei like 
Sn and Pb isotopes. Ca isotopes are expected to be inter- 
esting cases where the low-energy modes could eventually 
start to be more collective (with respect to light nuclei) 
and the nature of the low-lying excitations in terms of 
collectivity and fragmentation may be investigated. 



II. BRIEF SUMMARY OF THE FORMAL 
ASPECTS OF SRPA 

The SRPA equations are known since many years 
and have been derived by following different procedures, 
such as the equations-of- motion method [33|, the small- 
amplitude limit of the time-dependent density matrix 
method [40 . 41] , and a variational procedure introduced 
by da Providencia (42|. The main properties are also 
recalled in more recent works [3a, [S^l ■ 

The excited states in SRPA are superpositions of 1 
particle-1 hole [Iplh) and 2p2h configurations. The 
SRPA equations can be written in the compact form 



A 
-B* 



B 

-A* 



where: 



A 



X" 



Axi Ai2 
A21 A22 



xr 

x^ 



X" 



,B 



X" 



Bii B12 
B21 B22 



(1) 



y 



Y2^ 



In the above equations, '1' and '2' stand for Iplh and 
2p2h, respectively. Thus, An and Bn represent the 
usual RPA matrices, whereas the matrices A12 and B12 
couple Iplh with 2p2h configurations and the matrices 
A22 and 1^22 couple among themselves 2p2h configura- 
tions. The detailed expressions of these matrices can be 
found for example in Ref. [37[. If a density-dependent 
interaction like the Skyrme force is employed, rearrange- 
ment terms appear in the residual interaction. The usual 
RPA rearrangement terms appear in the matrices An 
and Bii. New types of rearrangement terms have been 
obtained for the other matrix elements in Ref. i38| within 
a variational derivation of the SRPA equations. The ex- 
pressions of these rearrangement terms are reported in 
Ref. f3S| and are used in this work. 

Before analyzing the results, we recall here two 
main properties of SRPA: the quasiboson approximation 
(QBA) is adopted in standard SRPA and is adopted here; 
the energy-weighted sum rules (EWSRs) are satisfied in 
SRPA as demonstrated formally in Ref. [33] and verified 
numerically in Ref. [37i [. 



III. RESULTS FOR 



'CA 



The paper is organized as follows. In Sec. II the main 
formal aspects of the SRPA model are briefly recalled. 
In Sec. Ill the low- lying strength distributions are ana- 
lyzed for the nuclei ^°'^^Ca and the transition densities 
associated to some states are displayed. In Sec. IV some 
comments about the spurious state are presented. We 
draw our conclusions in Sev. V. 



Both Ca isotopes are stable but a neutron skin has 
been measured experimentally in '*^Ca by proton and 
electron scattering experiments [4j|. The proton (neu- 
tron) radii found within the SGII-Hartree-Fock model 
are equal to 3.37 (3.32) fm and 3.41 (3.55) fm for the 
nuclei ^"Ca and ''^Ca, respectively. The experimental 
low-lying dipole response has been recently analyzed in 
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FIG. 1: RPA dipole strength distribution for ''^Ca. 



the two isotopes [Sl and the development of a low-energy 
strength, between 5 and 10 MeV, has been observed in 
^*Ca. From the theoretical point of view, it has been 
found that relativistic and non-relativistic (Q)RPA mod- 
els are not able to well describe the low-lying response in 
^®Ca because they either do not provide the good exci- 
tation energies (too high energies) or do not predict the 
experimental fragmentation of the peaks. For example, in 
recent calculations performed with the relativistic RPA 
model no strength has been found in the response below 
the excitation energy of 10 MeV [4J]. The same kind 
of results is obtained in Skyrme SGII-RPA calculations 
(Fig. [T]). On the other hand, a reasonable agreement (en- 
ergies and fragmentation) with the experimental results 
has been found within the ETFFS model [y, [2l| where a 
quasiparticle-phonon coupling is included. 

We perform SRPA calculations in spherical symmetry 
for the two isotopes '*°Ca and ^®Ca with the Skyrme in- 
teraction SGII. The technical details of these calculations 
are reported in Ref. [33|- Differently from Ref. [33|, the 
full rearrangement terms [3g| are used here in the SRPA 
matrices. Because of the zero-range of the Skyrme inter- 
action, a natural energy cutoff is not provided. Different 
procedures to treat this problem may be envisaged for 
future studies (see, for instance, the exploratory work 
presented in Ref. [l^)- In this work, we have introduced 
an energy cutoff (ECUT) on the 2p2h configurations. 
By varying it from 40 to 60 MeV we have verified that 
a reasonable stability of the results is achieved around a 
cutoff of 50-55 MeV. The total B{E1) and EWSRs val- 
ues, integrated up to an energy of 10 MeV, are shown in 
Table U for the isotopes ''"Ca and '**Ca as a function of 
the energy cutoff ECUT on the 2p2h configurations. 

The B{E1) distributions for different choices of the 
energy cutoff are plotted in Fig. [5] up to an excitation 
energy of 10 MeV for the nucleus ''^Ca. The employed 
transition operator is 
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FIG. 2: (Color online) SRPA dipole strength distribution for 
the ''^Ca isotope for increasing values of the energy cutoff 
ECUT (MeV) on the 2p2h configurations included in the cal- 
culations. 





^«Ca 


^"Ca 


ECUT 


EB(si) 


EWSRs 


Y:,B{E1) EWSRs 


40 


0.184 


1.623 


0.009 0.091 


45 


0.218 


1.895 


0.002 0.022 


50 


0.226 


1.944 


0.015 0.139 


55 


0.240 


2.049 


0.025 0.237 


60 


0.230 


1.964 


0.023 0.211 



TABLE I: Total B{El) {e^fm^) and EWSRs {e^fm^MeV) 
integrated up to 10 MeV as a function of the energy cutoff 
ECUT {MeV) on the 2p2h configurations for ''^Ca and ''°Ca. 
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where Cp and e„ are the kinematic charges, Cp = Ne/A 
and e„ = Ze/A, respectively. One can observe that the 
results do not change strongly starting from a cutoff of 
45 MeV. In what follows, we will analyze the results ob- 
tained with a cutoff of 60 MeV (bottom panel of Fig. ^ . 
These results are qualitatively of the same type as those 
found in Ref. p with the ETFFS approach which is also 
a beyond-mean-field model where the coupling is done 
with collective phonons instead of 2p2h configurations 
(as is done in SRPA). We can compare the location of 
the theoretical peaks with the experimental distribution 
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FIG. 3: (Color online) As in Fig. [2] but for the '^''Ca isotope. 
Please note that a different scale has been used in the ordinate 
with respect to Fig. [21 



(see, for instance, Fig. 2 of Ref. [g). We can distin- 
guish two regions: from 6 to 8 MeV and from 8 to 10 
MeV. Experimentally, the highest peak in the first re- 
gion is found at ~ 7 MeV whereas in our case we have 
several small peaks between 6 and 8 MeV and the highest 
peak is located around 6.2 MeV. In the interval between 
8 and 10 MeV, our response is more fragmented than 
the experimental one. Experimentally, peaks are found 
around 8.5, 9 and 9.5 MeV. Our highest peak is located 
at ^ 9.1 MeV. Finally, we have evaluated the response in 
the nucleus '"^Ca. Experimentally, a negligible strength 
has been found Jor this nucleus between 5 and 10 MeV 
(Fig. 2 of Ref. [6|). Our results for '^°Ca are displayed in 
Fig. [3] for cutoff values varying from 40 up to 60 MeV. 
Some peaks are actually found at low energy but the cor- 
responding strength is much lower than in ^^Ca (notice 
the different scales in the two figures). 

An interesting information that can be analyzed in con- 
nection with the strength distribution is the composition 
of the excitation modes in terms of Iplh and 2p2h con- 
figurations. By extracting the expression of A^i from the 
SRPA normalization condition. 
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we plot in Fig. |4]the B{E1) values corresponding to a 
cutoff of 60 MeV (upper panel, same as in bottom panel 
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FIG. 4: (Color online) For each state the B{E1) value corre- 
sponding to a cutoff of 60 MeV (upper panel) and the total 
Iplh contribution A'^i to the norm of the state defined in Eq. 
(|3]) (lower panel), are shown for **Ca. 



of Fig. [2|) and the quantity Ni (lower panel) for each 
excitation of the nucleus ^^Ca. One observes that all 
the excitations present a mixing of Iplh and 2p2h con- 
figurations. Those which have the highest Iplh content 
(around 50%) may be interpreted as excitations that al- 
ready exist in the RPA spectrum at higher energies (the 
first excitations in RPA are located around 11 MeV) and 
that are shifted down to lower energies due to the cou- 
pling with 2p2h configurations. In the SRPA spectrum 
we also see several states which present a dominant 2p2h 
nature and show a relatively large B{E1) despite their 
very low content of Iplh configurations (see for instance 
the energy region around 9 MeV). 

As already mentioned in the Sec. I, low-energy excita- 
tions in light nuclei are mostly single-particle excitations 
and cannot be interpreted as a collective motion of a 
skin against a core. In Ca isotopes, which are interme- 
diate cases between light and heavy nuclei, the nature of 
these excitations has not yet been clearly elucidated. In 
order to have a deeper insight into the properties of these 
low-energy modes we consider in more detail the states 
located at 5.95, 6.19 and 6.60 MeV as well as those lo- 
cated at 9.09 and 9.23 MeV. We can see that the first 
three states have a quite large and similar Iplh compo- 
nent but show a very different B{E1) value. In particu- 
lar, one notices that the state located at 5.95 MeV has 
almost no strength. By comparing among themselves the 
states with energies 6.19 and 6.60 MeV, one observes that 
Ni is larger while the B{E1) value is about one half in 
the second state with respect to the first one. A different 
description is provided for the states lying at 9.09 and 
9.23 MeV which are both almost entirely composed by 
2p2h configurations. In spite of the fact that the Iplh 
content is higher in the second state with respect to the 
first one, we notice that the second state has almost no 



E (MeV) 


B{E1,T ^0) 


B{E1,T = 1) 


5.95 


0.002 


0.013 


6.19 


0.176 


0.157 


6.60 


0.011 


0.036 


9.09 


0.099 


0.186 


9.23 


0.042 


0.012 



Neutron 

Proton 



Isoscalar 

Isovector 



TABLE II: Isoscalar B{E1,T = 0) {e^ frn^) and isovector 
B{E1,T = 1) {e^fm^) transition probabilities obtained by 
integrating the corresponding curves shown in Figs [5] and [6] 
multiplied by r. 



strength while the first state is the most collective in the 
low-lying spectrum. 

The nature of the low-energy peaks can be better ana- 
lyzed by looking at the associated transition densities. In 
Fig. [S] we compare the transition densities correspond- 
ing to the peaks of 5.95 (upper panels), 6.19 (middle 
panels) and 6.60 (lower panels) MeV. For each state, 
we show separately the neutron 6pn{r) and the proton 
Spp{r) (left), the isoscalar Spn{r) + Spp{r) and the isovec- 
tor Spn{r)~Spp{r) (right) transitions densities. The same 
quantities are shown in Fig. IH] for the states located at 
9.09 (upper panels) and 9.23 (lower panels) MeV. 

In Table |lT] we report for each state the isoscalar 
B{El,T = 0) and isovector B{El,T = 1) transition 
probabilities obtained by integrating the corresponding 
sets of Figs [5] and [5] multiplied by r. By looking at the 
transition densities of the first three states (Fig. [5]) one 
does not see any clear signature of an oscillation of the 
neutron skin against the core at the surface of the nu- 
cleus. On the contrary, especially in the external part of 
the nucleus, the protons and neutrons oscillate in phase. 
For the state located at 5.95 MeV, the isoscalar and 
isovector transition densities strongly oscillate giving al- 
most vanishing B{E1,T = 0) and B{E1,T = 1) values 
(see Table HI)) . A similar behavior is found for the state 
of 6.60 MeV. For the state lying at 6.19 MeV the cancel- 
lations are less important. Strong cancellations occur for 
the state located at 9.23 MeV, (lower panel of Fig. [SJ, 
resulting in very small B{E1,T = 0) and B{E1,T = 1) 
values (Table [TTI . A different situation is found for the 
most collective state located at 9.09 MeV. We see that 
the neutron transition density dominates over the proton 
one that is almost vanishing in the external part of the 
nucleus while in the interior the two densities oscillate 
out of phase. This leads to a strong mixing of isoscalar 
and isovector components. 

Another interesting analysis that can be done for these 
excitation modes is related to their collectivity in terms 
of number and coherence of the different Iplh configu- 
rations which contribute the total transition probability. 
We present an analysis similar to that done in Ref. |17| . 
In SRPA as well as in RPA the reduced transition prob- 
ability for a one-body operator describing the excitation 
from the ground state to a state v can be written as 
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FIG. 5: (Color online) Neutron and proton transition densi- 
ties (left) and the corresponding isoscalar and isovector ones 
(right) associated to the peaks located at 5.95 (upper panels), 
6.19 (middle panels) and 6.60 (lower panels) MeV. 





E-- 


= 5.95 MeV 






ph conf. 


E (MeV) 


Aph 


bph{El) 


Fp\ 


(2p3/2, 1^5/2)'^ 


17.499 


0.001 


0.021 


1.343 


(1/7/2, 1^5/2)" 


11.732 


0.062 


-0.470 


3.304 


(2^3/2,281/2)" 


12.444 


0.022 


0.140 


1.726 


(2pi/2,2Si/2)" 


14.133 


0.004 


0.043 


-1.233 


(2p3/2,ld3/2)" 


12.120 


0.172 


0.107 


0.451 


(2^1/2,1^3/2)" 


13.809 


0.007 


-0.050 


1.053 


(1/5/2, ld3/2)" 


13.867 


0.012 


0.188 


2.756 


(2^3/2,281/2)" 


11.773 


0.006 


0.060 


1.737 


(2^1/2,281/2)" 


13.556 


0.001 


-0.011 


-1.248 


(2p3/2, 1^3/2)" 


10.329 


0.145 


-0.071 


0.456 


(2pi/2, 1^3/2)'^ 


12.112 


0.001 


0.014 


1.084 


Partial Sum 




0.433 


-0.031 




Total Sum 




0.435 


-0.054 





TABLE III: Particle-hole configurations which give the ma- 
jor contributions to the dipole low-lying state located at 5.95 
MeV. For each ph configuration, the energy, the contribution 
to the norm of the state Aph, the partial contribution to the 
reduced transition amplitude hph (e fm) and the matrix ele- 
ment of the transition operator Fp^ are reported. The super- 
scripts TT, V refer to proton and neutron states, respectively. 
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FIG. 6: (Color online) As in Fig. [5] but for the states located 
at 9.09 and 9.23 MeV (upper and lower panels, respectively). 





E-- 


= 6.19 MeV 






ph conf. 


E (MeV) 


Aph 


bpH{El) 


FX 


(1/7/2, 1^5/2)" 


11.732 


0.156 


-0.727 


3.304 


(2p3/2,2Sl/2)'' 


12.444 


0.148 


0.373 


1.726 


(2pi/2,2Si/2)" 


14.133 


0.003 


-0.043 


-1.233 


(2^3/2,1^3/2)" 


12.120 


0.003 


-0.015 


0.451 


(2^1/2,1^3/2)'" 


13.809 


0.073 


0.166 


1.053 


(1/5/2,1^3/2)" 


13.867 


0.018 


0.250 


2.756 


(2p3/2, 1^5/2)'' 


15.683 


0.000 


0.012 


1.410 


(2^3/2,281/2)" 


11.773 


0.015 


-0.079 


1.737 


(2p3/2, 1^3/2)'' 


10.329 


0.040 


0.038 


0.456 


(2pi/2, 1^3/2)" 


12.112 


0.001 


-0.016 


1.084 


(159/2,1/7/2)" 


11.364 


0.003 


-0.106 


4.171 


Partial Sum 




0.461 


-0.147 




Total Sum 




0.465 


-0.163 





TABLE IV: As in Table HII] but for the state located at 6.19 
MeV 





E-- 


= 6.60 MeV 






ph conf. 


E (MeV) 


Aph 


hpH{El) 


FX 


(2p3/2, Ids/a)" 


17.499 


0.001 


0.031 


1.343 


(1/7/2,1^5/2)" 


11.732 


0.012 


0.226 


3.304 


(1/5/2,14/2)" 


19.246 


0.001 


0.010 


-0.734 


(2p3/2,2Si/2)" 


12.444 


0.190 


0.421 


1.726 


(2pi/2,2Si/2)" 


14.133 


0.002 


-0.049 


-1.233 


(2^1/2,1^3/2)" 


13.809 


0.335 


-0.350 


1.053 


(1/5/2, 1^3/2)" 


13.867 


0.014 


-0.193 


2.756 


(2^3/2,281/2)" 


11.773 


0.005 


-0.044 


1.737 


(2^3/2,1^3/2)" 


10.329 


0.035 


0.035 


0.456 


(2^1/2,1^3/2)" 


12.112 


0.018 


-0.060 


1.084 


(2^5/2, 1/7/2)" 


12.698 


0.001 


0.015 


1.212 


(3^5/2,1/7/2)" 


15.625 


0.000 


0.010 


1.055 


(1^9/2,1/7/2)" 


11.364 


0.001 


0.058 


4.171 


Partial Sum 




0.615 


0.110 




Total Sum 




0.622 


0.104 





TABLE V: As in Table [XITl but for the state located at 6.60 
MeV 



where i^A are the multipole transition amplitudes as- 
sociated to a Iplh configuration. We remark that also 
in the case of SRPA only Iplh amplitudes appear in the 
expression of the transition probability. A different situa- 
tion would occur if a two-body operator is considered. In 
Refs. [231 and [33, for example, the study of the double 
giant dipole resonance in ""^Ca and ^^O has been carried 
out by using a two-body operator. In the spirit of a mul- 
tiphonon picture, the latter is built as a product of two 
one-body dipole operators. As shown in Fig. |4l the low- 
lying dipole states that we obtain in SRPA have a strong 
2p2h nature. It could thus be interesting to investigate 
their properties by using a two-body transition operator. 
In particular, the use of a transition operator contain- 
ing both one-body and two-body terms is expected to 
affect the strength distribution and eventually the total 
strength associated to this energy region. On the other 
hand, it is not clear which kind of one-body multipole 
operators should be taken into account here to construct 
the two-body operator. This investigation is left as a 
subject for a future work. 



In Tables ITTTll VIII we report the particle- hole configura- 
tions which provide the major contributions to the dipole 
modes for the five states analyzed in Figs [5] and [51 For 
each configuration we report the unperturbed energy, the 
contribution Api^ to the norm of the state. 



Anh — A„i, 



ph\ 



\^ph\ , 



with 



ph 



(5) 



the partial contribution hph to the reduced transition am- 
plitude (see Eq. ^) and the matrix element of the tran- 
sition operator. 

We first briefly discuss the case of the RPA IVGDR 
whose collective features are well known. In Fig. [7] the 
partial contributions hph corresponding to each \p\h con- 
figuration for the first peak located at 17.33 MeV (Fig. 
[T]) are shown. The bars correspond to each value of bph 
associated to a single configuration while the continu- 
ous line is the cumulative sum of the contributions. The 
dashed line separates the proton from the neutron con- 
figurations which are ordered according to their increas- 



ing energy. We can clearly see that the contributions of 
many proton and neutron Iplh configurations sum up 
coherently to provide the total B{E1). This coherent be- 
havior of protons and neutrons is due to the minus sign 
in the definition of the isovcctor transition operator, Eq. 
^ , and is not in contrast with the isovcctor character of 
this excitation where neutrons and protons oscillate out 
of phase. 

In Figs. [8] and [9] the partial contributions bph corre- 
sponding to each Iplh configuration for the excitation 
modes located at 5.95, 6.19 and 6.60 MeV (Fig. El and 
at 9.09 and 9.23 MeV (Fig. 9) are plotted. 

In the upper panel of Fig. [8] we see that, for the state 
lying at 5.95 MeV, several Iplh configurations present 
non negligible bph values. This is true especially for the 
proton configurations and is also indicated by the behav- 
ior of the proton transition density in Fig. [5] However, 
the total amplitude is very small since strong cancella- 
tions occur. The same holds for the state located at 6.19 
MeV (middle panel of the same figure), but only for the 
proton Iplh configurations whose cumulative sum is al- 
most zero. The total transition amplitude is given in 
this case only by the neutron configurations. A differ- 
ent result is found for the third state shown in the lower 
panel, where the strong cancellation occurs for the neu- 
tron configurations. For the highest state lying at 9.23 
MeV (lower panel of Fig. [5]) we observe strong cancel- 
lations for both neutrons and protons leading to a very 
small total transition amplitude. A more interesting sit- 
uation is obtained for the most collective state, located 
at 9.09 MeV (upper panel of Fig. [5]). We observe also 
in this case a strong cancellation of the proton contribu- 
tions while a quite coherent behavior is exhibited by the 
neutron Iplh configurations. In particular, we observe 
strong contributions coming from the outermost neutrons 
(see also Table FVl)) . This result, together with the profile 
of the corresponding transition density indicates that this 
state shows some features usually associated to pygmy 
resonances. 

In TableEIIIlwe compare the total B{E1) and EWSRs 
integrated up to 10 MeV and the corresponding centroid 
energies with the experimental values ^] for the nuclei 
4°Ca and "^Ca. We see that the SRPA total B{E1) is 
much larger, almost by a factor 4, than the experimental 
value. This is due to the larger number of states obtained 
in SRPA with respect to the experimental spectrum and, 
at the same time, to their higher strength (Fig. [T0| . 
The same kind of discrepancies is found for the EWSRs 
while the SRPA centroid energy is very close to the ex- 
perimental value. Regarding these large differences some 
comments are in order. First we recall that, because of 
the non-local terms of the Skyrme interaction, the dou- 
ble commutator sum rule is enhanced with respect to the 
classical sum rule by a factor 1.35 for SGII. However, 
this is not enough to explain the strong deviation with 
respect to the experimental values. It could be also inter- 
esting to analyze whether or not this kind of discrepancy 
may depend on the choice of the Skyrme interaction. As 
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FIG. 7: (Color online) Partial contributions bph of the re- 
duced transition probability versus the ordering number of 
the Iplh configurations for the first peak of the RPA IVDGR 
located at 17.33 MeV (Fig. [1]). The dashed line separates the 
proton from the neutron configurations. The configurations 
are ordered according to their increasing energy. The bars 
corresponds to the individual bph contributions while the full 
red line is the cumulative sum of the contributions. 
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FIG. 8; (Color online) As in Fig. [7| but for the state at at 
5.95 MeV, 6.19 MeV and 6.60 MeV. 





E-- 


= 9.09 MeV 






ph conf. 


E (MeV) 


Aph 


bMEl) 


fX 


(2^3/2,1^5/2)'" 


17.499 


0.002 


-0.039 


1.343 


(1/7/2,14/2)" 


11.732 


0.001 


-0.098 


3.304 


(1/5/2,1^5/2)" 


19.246 


0.006 


-0.032 


-0.734 


(2pi/2,2Si/2)- 


14.133 


0.005 


0.045 


-1.233 


(2p3/2, 1^3/2)" 


12.120 


0.007 


-0.023 


0.451 


(2pi/2,ld3/2)" 


13.809 


0.000 


-0.014 


1.053 


(1/5/2,1^3/2)" 


13.867 


0.017 


0.197 


2.756 


(2p3/2,ld5/2)" 


15.683 


0.006 


-0.044 


1.410 


(2^3/2,281/2)" 


11.773 


0.006 


-0.055 


1.737 


(2pi/2,2Si/2)" 


13.556 


0.001 


-0.015 


-1.248 


(2^3/2,1^3/2)" 


10.329 


0.011 


-0.020 


0.456 


(2pi/2, 1^3/2)'' 


12.112 


0.001 


0.018 


1.084 


(1/5/2,1^3/2)" 


14.072 


0.009 


0.105 


2.675 


(2^5/2,1/7/2)" 


12.698 


0.001 


0.014 


1.212 


(I59/2, l,/7/2)" 


11.364 


0.008 


0.149 


4.171 


Partial Sum 




0.082 


0.189 




Total Sum 




0.083 


0.190 





TABLE VI: As in Table [Till but for the state of energy 9.09 
MeV 





E-- 


= 9.23 MeV 






ph conf. 


E (MeV) 


Aph 


bph{El) 


Fpi 


(2p3/2,ld5/2)" 


17.499 


0.002 


0.033 


1.343 


(1/7/2,1^5/2)" 


11.732 


0.001 


-0.059 


3.304 


(2p3/2,2Sl/2)" 


12.444 


0.000 


-0.012 


1.726 


(2pi/2,2Sl/2)" 


14.133 


0.001 


0.027 


-1.233 


(2p3/2, 1^3/2)" 


12.120 


0.023 


0.040 


0.451 


(2p3/2,ld5/2)" 


15.683 


0.008 


-0.055 


1.410 


(2p3/2,2Si/2)" 


11.773 


0.103 


-0.226 


1.737 


(2pi/2,2Si/2)" 


13.556 


0.001 


0.019 


-1.248 


(2p3/2, 1^3/2)" 


10.329 


0.004 


0.012 


0.456 


(2pi/2,ld3/2)" 


12.112 


0.056 


0.106 


1.084 


(2d5/2, 1/7/2)"^ 


12.698 


0.001 


0.013 


1.212 


(lff9/2,l,/7/2)" 


11.364 


0.001 


0.052 


4.171 


Partial Sum 




0.201 


-0.050 




Total Sum 




0.202 


-0.037 





TABLE VIL 
MeV 



As in Table IIIII but for the state lying at 9.23 



a check, we have performed SRPA calculations by using 
the parametrization SLy4 and the same kind of devia- 
tions has been found. In a recent work [4g| it has been 
suggested that the low-lying dipole strength distribution 
could be be eventuahy related to the slope of the symme- 
try energy. This kind of analysis is however beyond our 
present scopes and it will be performed in future investi- 
gations. Finally, we recall that theoretical B{E1) values 
much larger than the corresponding experimental results 
have also been found within the ETFFS model for the 
nucleus '*'*Ca; it has been shown that this discrepancy is 
related to the use of some approximations in the treat- 
ment of pairing and continuum coupling [21| . We also 
mention that in previous experimental measurements a 
much larger strength had been found in the energy re- 
gion from 5 to 10 MeV [43|. Finally, as mentioned above, 
since the low-lying dipole states obtained in SRPA have 
a strong 2p2h nature, it could more appropriate the use 
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FIG. 9: (Color onhne) As in Fig. [8] but for the state at 9.09 
and 9.23 MeV (upper and lower panel, respectively). Please 
note that a different scale has been used in the ordinate with 
respect to Fig. [HI 







^«Ca 


^"Ca 


j:b{ei) 

(10-» e^ fm^) 


SRPA 
Exp 


230 
68.7 ± 7.5 


23 
5.1 ± 0.8 


J:^EA{E1) 
10"^ e^ fm^ MeV 


SRPA 
Exp 


1964 
570 ± 62 


211 
35 ± 5 


J-'centroid 

MeV 


SRPA 
Exp 


8.54 
8.40 


9.17 
6.80 



TABLE VIII: Total B{E1) and EWSRs integrated up to 10 
MeV and corresponding centroid energies obtained in SRPA 
compared with the experimental values |Q| for the ^'''^^Ca 
isotopes. 



of a transition operator containing both one-body and 
two-body terms. Of course, the use of such a more gen- 
eral operator would affect the total strength associated 
to this energy region. This investigation is however left 
as a subject for a future work. 



IV. MIXING WITH THE SPURIOUS STATE 

Some comments about the spurious state and its pos- 
sible mixing with the physical dipole modes are in order. 
The Thouless theorem on the EWSRs [4^ is very impor- 
tant in the framework of RPA and it holds also in SRPA 
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FIG. 10: Comparison of the experimental B(E1) strength dis- 
tribution [^ (upper panel) with the SRPA calculations (lower 
panel) for "^^Ca. 
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FIG. 11: (Color online) Dipole strength distribution for the 
**Ca isotope obtained in SRPA without (full red lines) and 
with a renormalizing factor Vmn = 0.91 in the residual inter- 
action (dashed blue lines). Both calculations are done with 
an energy cutoff on the 2p2h configurations of 60 MeV. 



[39|. It guarantees that spurious excitations correspond- 
ing to some symmetries separate out and are orthogonal 
to the physical states. This separation is obtained only in 
completely self-consistent calculations, that is, when the 
same interaction is used at both the HF and the RPA 
(or SRPA) level. In the case of dipole excitations, the 
center-of-mass motion should appear at zero energy and 
the EWSRs should be satisfied. Since the Coulomb and 
spin-orbit terms are not taken into account in the residual 
interaction, our calculations are not fully self-consistent 
and violations of the EWSRs are found (not larger than 
2 — 3%). A currently adopted procedure to estimate the 
mixing with spurious components consists in using the 
isoscalar one-body operator corrected for the center-of- 
mass motion. This has been done, for example, in the 
study of giant resonances in Ref. [33] ■ However, this 
check which is generally employed in RPA calculations 
is suitable to analyze excitations which are mainly com- 
posed by Iplh states. In our case, as shown in Fig. |3J 
Iplh and 2p2h components are strongly mixed and we 
thus choose a more appropriate procedure to estimate 
the mixing with spurious components. 

The residual interaction has been multiplied by a 
renormalizing factor to shift to zero the energy of the 
spurious mode. In our RPA calculations, the spurious 
state lies at about 3.5 MeV exhausting more than 95% 
of the isoscalar EWSRs; by using a renormalizing fac- 
tor of 1.09 its energy goes down to 0.2 MeV while the 
rest of the isoscalar and isovector distributions remains 
practically unaffected. In SRPA, as a consequence of 
the coupling with the 2p2h configurations, the spurious 
state is pushed down to lower energies with respect to 
RPA. Depending on the energy cutoff on the 2p2h con- 
figurations, its energy can be negative or imaginary in 
some cases. We stress that the use of the Hartree-Fock 
ground state that minimizes the energy guarantees that 



the RPA matrix is positive definite. This does not hold 
in SRPA. For the SRPA calculations with ECUT = 60 
MeV we have verified that, by using a renormalizing fac- 
tor of 0.91, the spurious mode is located at about 0.12 
MeV. The total B{E1) (EWSRs) integrated up to 10 
MeV does not change much when the renormalizing fac- 
tor is introduced: from 0.230 e^ fm^ (1.964 e^ fm^ MeV) 
to 0.221 e^ fm^ (1.905 e^ fm^ MeV) . In Fig. [TT] we 
show the dipole strength distribution obtained for the 
■**Ca isotope without (full red lines) and with (dashed 
blue lines) the inclusion of the renormalizing factor in 
the residual interaction. Both calculations are done with 
an energy cutoff on the 2p2h configurations of 60 MeV. 
We see that the strength distribution is not strongly af- 
fected, the main difference being a shift of few hundreds 
KeV of some peaks and a change of their corresponding 
transition probability. This analysis seems to indicate 
that the presence of possible spurious components does 
not affect much our results. 



V. CONCLUSIONS 

In this work, we have analyzed the low-energy dipole 
spectrum (from 5 to 10 MeV) for the stable nuclei 
^°Ca and ^^Ca in the framework of the Skyrme-SRPA 
model. The Skyrme interaction SGII is used. Al- 
most no strength is found in this energy region for the 
nucleus *°Ca whereas a non negligible strength is ob- 
tained for the neutron-rich nucleus '^^Ca. The distribu- 
tion and the fragmentation of the peaks is in reasonable 
agreement with the corresponding experimental measure- 
ments. This kind of results cannot be provided by the 
standard RPA model: SGII-RPA calculations do not lead 
to any strength in the energy region from 5 to 10 MeV. 
However, we have found a B(E1) value integrated up 



10 



to 10 MeV which is quite larger than the corresponding 
experimental result. 

The inclusion of 2p2h configurations and their coupling 
with the Iplh ones and among themselves within SRPA 
has a two-fold effect on the low-lying dipole strength: (i) 
states that already exist in RPA are shifted to lower en- 
ergies. These states maintain a quite strong Iplh nature; 
(ii) several other states of almost pure 2p2h character ap- 
pear. They are excited by the one-body dipole operator 
through their Iplh components. In the nucleus "^^Ca one 
of these states shows the largest transition strength and 
displays some features generally associated with a dipole 
pygmy resonance. 

A detailed analysis of the main excitations that com- 
pose the strength distributions is done: the content of 
Iplh and 2p2h configurations is studied for some peaks. 
The transition densities are shown and the collectivity of 
the peaks is investigated in terms of the number and the 
coherence of the single configurations that mainly con- 



tribute. As far as the transition densities are concerned, 
one may conclude that in general they do not display the 
typical profile well known for pygmy resonances, except 
for the state located at 9.09 MeV. We have also observed 
that about 10-12 configurations contribute to each peak. 
In this sense one can say that there is some collectivity, 
however, strong cancellations occur in most cases and the 
single configurations do not sum up in a coherent way. 

Our main conclusion is that, even if a low- lying dipole 
response is found experimentally and predicted theoreti- 
cally in the isotope ^^Ca, one cannot really describe these 
excitations as pygmy resonances except for the most col- 
lective peak located at 9.09 MeV. This suggests that this 
nucleus is still too light to present clear signatures of an 
oscillation of the neutron skin against the internal core 
and that individual degrees of freedom are still dominant 
in the description of the dipole low-energy spectrum as 
for lighter nuclei. 
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